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Shock waves in cold nuclear matter, e.g., those induced by a collision of two neutron stars, can 
generate a large number of gamma photons via parametric resonance. We study the resonant 
production of gamma rays inside a shocked neutron star and discuss the possible astrophysical 
consequences of this phenomenon. 
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I. INTRODUCTION 

Neutron star collisions rank among the most energetic 
events expected to take place in the Universe, which 
makes them a natural candidate for the source of ob- 
served gamma-ray bursts (GRB) 0. The discovery of 
the afterglow j!| associated with some of the GRB and 
the isotropy of the GRB both support the hypothesis that 
the GRB originate at cosmological distances with a red- 
shift of order one. The total energy ~ 10 53 erg released 
when the two neutron stars merge is sufficient to explain 
the GRB, although it is not clear what fraction of that 
energy is converted to gamma rays. Another puzzling 
feature of the GRB is their non-thermal spectrum. 

Understanding the physics of the explosion that fol- 
lows the collision of two neutron stars is of great im- 
portance and is the subject of intense studies The 
approach and the early stages of the interaction of the 
two neutron stars are accompanied by powerful acous- 
tic shock waves that propagate through nuclear matter 
and eventually dissipate their energy into heat. During 
this period of time, the repetitive superconducting phase 
transitions can take place in part of the star's volume due 
to the sharp dependence of the proton energy gap on den- 
sity. The relaxation of the proton condensate to the po- 
tential minimum can be accompanied by a non-thermal 
resonant production of gamma rays Q in the MeV en- 
ergy range. If the parametric resonance is efficient, the 
power transferred to the coherent gamma quanta may ex- 
ceed 10 52 erg/s. Eventually, the nuclear matter is heated 
above the critical temperature ~ 0.5 MeV and the pro- 
duction of the gamma rays comes to a halt. 

The photons produced inside a neutron star cannot de- 
cay through pair production 7 — > e + e~ because the de- 
generate electrons have the chemical potential in excess of 
100 MeV, far greater than the photon energy. The decay 
into the electron-positron pairs is, therefore, prohibited 
by the Pauli exclusion principle. The photons undergo a 
Compton scattering off the electrons near the Fermi sur- 
face which is sufficiently strong to keep the photons from 



escaping. While such scattering (or comptonization) may 
change the spectrum of the gamma-component, the num- 
ber of photons remains the same. Since the extreme 
electron degeneracy is maintained even at the highest 
temperatures achieved in the fireball ||, the gamma-ray 
component is present in the nuclear matter at the time 
the latter is dispersed by the explosion. 

The latter can have several consequences. In partic- 
ular, 7-quanta, abundantly present in nuclear matter at 
the onset of the explosion, can be released when the fire- 
ball erupts. The detailed investigation of this signal lies 
outside the scope of this paper. However, it is plausi- 
ble that the gamma-rays emitted at that point would 
have a non-thermal spectrum. At later times, when 
the fireball reaches a high temperature, other sources of 
gamma emission become dominant. We predict, there- 
fore, a qualitative difference in the spectrum of gamma- 
rays emitted during the first milliseconds of the collision. 

In this paper we will not attempt to understand the 
emission of gamma-rays from the fireball. Instead, we 
will concentrate on the phenomenon of resonant produc- 
tion of photons which transfers a fraction of the gravita- 
tional energy into the gamma-quanta inside the neutron 
star. This process is of fundamental interest on its own 
and it may have important consequences. 



II. BASIC IDEA 

Nuclear matter in the interior of a neutron star may 
be superconducting. The existence of the superfiuid pro- 
ton condensate depends on several theoretical assump- 
tions, some of which may be hard to justify. The proper 
usage of the many-body techniques and the choice of 
macroscopic degrees of freedom are by no means obvi- 
ous. However, assuming that the theoretical framework 
ofRefs. §-§ is valid, one can elucidate some generic fea- 
tures of proton superfluidity. The most important one 
for us is that the energy gap depends sharply on den- 
sity, as shown in Fig. |]. An acoustic shock wave passing 
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FIG. 1. The energy gap as a function of matter density 
inside a neutron star j£j. Superconductivity occurs in some 
range of densities, I0 14 < p < 6 x I0 14 g/cm 3 . 



through the star can produce significant changes in the 
density and cause repetitive superconducting phase tran- 
sitions occurring with periods of order the acoustical time 
scale T a ~ I0~ 7 to 10~ 3 s Q. The relaxation time of the 
proton condensate is of order MeV -1 , or 10~ 20 s <C r tt , 
and therefore the system quickly settles in the potential 
minimum after the passage of each shock wave. 

In the presence of a magnetic field the superconduct- 
ing phase transition is first-order |^] and occurs in two 
stages. First, a bubble of the superconducting phase nu- 
cleates and expands; then the proton condensate may 
oscillate around the minimum. If the magnetic field is 
close to critical, B ~ B c , the two minima in the potential 
shown in Fig. |2j are nearly degenerate and no subsequent 
coherent motion of the condensate <f> takes place after 
the transition through bubble nucleation. However, for 
a smaller magnetic field, the so called "escape point" <f> e 
is different from the global minimum 0o & n d the scalar 
condensate may oscillate around the minimum. 

During the oscillations of the order parameter a pho- 
ton has effectively a time-dependent mass proportional 
to the value of <f>. This may, in some cases, signal a co- 
pious production of gamma-quanta through parametric 
resonance |1 . In this paper we study the parametric res- 
onance both analytically and numerically for different for 
some sample values of nuclear density and magnetic field 
that ranges from 10 12 to 10 15 G. 



III. POTENTIAL AND EQUATION OF MOTION 

A natural framework for describing the relaxation of 
the proton condensate after the phase transition is time- 
dependent Ginzburg-Landau theory jT^Jll[]. The appli- 
cability of such description is limited to cases where the 
Joule heat loses are small E2]. Otherwise, the interac- 
tions become essentially non-local and a simple equation 
of motion for a scalar field seizes to be valid. 



Here we will assume the validity of the time-dependent 
Ginzburg-Landau theory |ll],[l2| and will apply it to the 
superconducting proton condensate inside a neutron star. 
One should take this approach with a grain of salt given 
the lack of empirical knowledge with respect to the pro- 
ton superconductivity in nuclear matter. However, we 
hope that - crude an approximation this may be - it 
may help identify the main features of the resonant pro- 
duction of gamma-rays. 

We write the equation of motion for the order param- 
eter as 
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where £f = pp/^m^ is the Fermi energy of the proton 
condensate and c = (28£(3)/3"7r 3 ) ep /T c is a constant 
characterizing the condensate. The friction term then 
has a magnitude 8e/3c ~ 7.37T C ~ 4.2Ao, of order a few 
MeV, which is comparable to the oscillating frequency u> 
of the condensate around its minimum </>o , and therefore 
cannot be ignored. On the other hand, the gradient term 
in Eq. (Q) is negligible in our case Q, and thus we are 
dealing with a (locally) homogeneous condensate. 

The effective potential for the order parameter 4> can 
bc determined from the physical properties of the conden- 
sate. After the acoustic wave has restored the symmetry, 
the energy density subsides and a new superconducting 
phase transition takes place. The difference in energy 
density associated with the proton condensate is p3| 
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where Ao is the proton energy gap, and m* is the effective 
proton mass in nuclear matter, which is somewhat lower 
than the bare proton mass m p . The equilibrium value 
of the order parameter is given by 4>q — n p /2m,, where 
n p = Ypp/rrip is the number density of protons inside the 
neutron star, which make up a fraction Y p 0.03 of all 
baryons. 

The way the superconducting phase transition occurs 
depends very strongly on the presence of a magnetic field 
B. Such a field creates a barrier in the effective potential 
between the symmetric and the superconducting phases. 
The height of the barrier is approximately given by the 
total energy density in the magnetic field, 
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This barrier makes the phase transition first-order, with 
creation of bubbles of superconducting phase. After the 
bubble nucleation, the proton condensate oscillates and 
soon settles (due to friction) at its equilibrium value 4> - 
We have shown in Fig. || the effective potential U(4>) 
for a region inside the neutron star with matter den- 
sity p = 10 14 g/cm 3 and magnetic field B = 10 15 G, 
which corresponds to A ~ 0.35 MeV, O ^ 50.3 MeV, 
AU((j)o) =s 3000 MeV 4 and A17(&) ~ 200 MeV 4 . 
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FIG. 2. Effective potential for the proton condensate, in 
units of MeV 4 . The small bump at <f>\ is due to a magnetic 
field B. 
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FIG. 4. The growth parameters /j,k as a function of mo- 
menta in units of the effective frequency, k/u>, after one oscil- 
lation of the condensate. 
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FIG. 3. The oscillations of the order parameter associ- 
ated with the proton condensate around the minimum of the 
effective potential. Time is given in units of the effective fre- 
quency. Due to the friction term, the amplitude of the proton 
condensate diminishes in a few oscillations. 



Depending on the relative size of the friction term, 
the condensate settles at <j>o after a few oscillations, or 
without oscillations at all. We have shown in Fig. || the 
oscillations of the order parameter (j> in the case of the 
parameters of Fig. 0, which give a friction coefficient of 
order 4.2 A ~ 0.31 u and w ~ 4.7 MeV. This allows a 
few oscillations to take place, which is crucial for the res- 
onant production of photons, as we will see in the next 
Section. 



IV. RESONANT PHOTOPRODUCTION 

The proton condensate comprises Cooper pairs of 
charge 2e, which couple to the electromagnetic field 
through the term (2e) 2 </> 2 A^A^ in the unitary gauge, in 
which the scalar field <f> is real. As the condensate os- 
cillates around cf>o, it induces an effective periodic mass, 



which in some cases may stimulate parametric resonant 
production of photons Q . The mechanism is similar to 
that discussed in connection with axion clumps p4| , as 
well as preheating after inflation Ea], where explosive 
production of bosons may occur under special circum- 
stances 

We will write the gauge field as A^(x) = x(x)e l _ l (x), 
where e M (cc) is a polarization vector and x( x ) can be ex- 
panded in Fourier modes, Xk(t), that satisfy the evolution 
equation 

X fc + (fc 2 + 2(2e)VW)xfc(*) = 0, (4) 
with an effective mass proportional to the condensate 

(j>(i) ~ 0o (l + $exp(-earf/2) sinwt) , (5) 



where $ = e /0o and e ~ 4.2Ao/w is the deca 
stant for the condensate oscillations. Equation 



be written as a Mathieu equation |17| with coefficients 
(z = ut/2) 



con- 
can 
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q{z) = 4gr $cxp(-ez) , 
q = 8e 



n — 327T Ct cm 



(6) 
(7) 

(8) 



Note that the effective parameter q{z) that determines 
the strength of the resonance decreases exponentially 
with time. If e is too large, the parametric resonance 
is weak, and this mechanism is inefficient. On the other 
hand, if e is small, then the condensate oscillates several 
times and causes an explosive production of gamma rays 
with energy of a few MeV. 

At density p w 10 14 g/cm 3 and magnetic field B ~ 
10 15 G, parametric amplification of the MeV photons can 
take place. In this case, e = 0.31 and Aq ~ 340. We 
have plotted pu as a function of k, after one oscillation 
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of the condensate, in Fig. |4|, from which we can deduce 
the spectrum 

n k ~ - exp(p k uit) . (9) 

Nuclear matter is opaque for photons with energy of or- 
der the electron chemical potential, which limits the spec- 
tral band in which the modes are amplified. The spec- 
trum n k has, therefore, an ultraviolet cut-off at the mo- 
mentum associated with the electron Fermi energy inside 
the neutron star, beyond which photons are absorbed and 
the resonance shuts down. This occurs at wavelengths 
k c ~ 100m e = 51 MeV. For the model discussed above, 
with frequency of oscillations uj = 4.7 MeV, this corre- 
sponds to k/ui ~ 11. 

Furthermore, after a few oscillations, backreaction oc- 
curs, that is, the number of created photons is so large 
that it dominates the frequency of oscillations of the pro- 
ton condensate, m? = u> 2 + 8e 2 (x 2 )- The backreaction 
sets in when [|15| 

«x 2 »=^ rw^^~cA do) 

2tH J uj k 8e 2 

where u 2 , — k 2 + 8e 2 <fi 2 (t) is the effective frequency of the 
mode k, and the integration is up to the physical cut-off, 
k c . For the model in hand, we find that backreaction 
takes place after about one oscillation, much before the 
friction term in Eq. (^) has significantly decreased the 
oscillation amplitude of the condensate. 

The total energy density in photons produced during 
the resonance is, therefore, 

1 f kc 

Ay(*) = / dkk 2 Lo k n k (t) , (11) 
27r Jo 

We have plotted this energy density in Fig. ||. p~/(t) 
reaches an asymptotic value after a few oscillations. Be- 
cause of backreaction, the energy density produced is 
that after the first oscillation, of order 10 5 MeV 4 . It 
may seem like a large value, but actually the energy den- 
sity produced via this mechanism is just a small fraction 
of the neutron star density, p 1 ~ 8 x 10 4 MeV 4 ~ 10~ 4 p. 

We have analyzed the resonant photoproduction nu- 
merically for magnetic field from 10 12 G to 10 16 G and 
found no significant deviation from the sample values de- 
scribed above. The main effect of the magnetic field is to 
produce a potential barrier at <fi = </>i, as shown in Fig. [|. 
A large magnetic field may, however, have some other 
effects on nuclear matter. For example, it can modify 
the particle composition, in particular the distribution 
of protons and electrons p8[ . This, in turn, can affect 
the resonant photoproduction. Other parameters, such 
as the electron chemical potential, the proton fraction 

, etc. ^ can differ significantly from the sample values 
we took. However, our numerical analyses show that a 
resonant production of gamma rays from the oscillating 
charged condensate is possible for a wide variety of pa- 
rameters. 
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FIG. 5. The energy density, in units of MeV 4 , of gamma 
ray photons produced via parametric resonance, as a function 
of time. 

V. IMPLICATIONS FOR THE OBSERVATION OF 
GAMMA-RAY BURSTS 

The MeV-energy gamma rays produced in each cycle 
of oscillations cannot decay the way they would decay in 
vacuum, through electron-positron pair production, be- 
cause the electrons are highly degenerate inside the neu- 
tron star. Pauli exclusion principle prevents a production 
of electrons with energies less than the electron chemi- 
cal potential, which is of order 100 MeV. Compton scat- 
tering off the electrons and protons near the Fermi sur- 
face is kinematicly suppressed but is not forbidden. The 
corresponding mean free path is A ~ 10~ 9 cm. Comp- 
tonization generally preserves the number of photons and 
tends to equalize the temperatures of electrons and pho- 
tons . Given enough time, the photons would diffuse 
out of the star. However, since the merger of the two 
neutron stars is characterized by very short time scales, 
all the gamma rays non-thermally produced by the co- 
herent oscillations of the proton condensate remain inside 
the star when the fireball erupts. At that point they can 
leak out and be observed as a short gamma-ray burst, or 
more likely, as a component of the GRB. 

The energy stored in the non-thermal bath of gamma 
ray photons inside the neutron star is Ej = p 7 V ~ 3 x 
10 50 erg. If the collision of neutron stars releases these 
photons within a time scale of order 1 ms < r < 0.1 s, 
then the power generated is of order 10 52±1 erg/s, which 
corresponds to the power emitted in the observed gamma 
ray bursts Q . 

The resonant photoproduction is active in a spherical 
shell with density 10 14 — 10 15 g/cm 3 that contains most 
of the neutron star mass (n]] , and where the acoustically 
driven superconducting phase transitions can take place. 
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VI. CONCLUSION 

Strong mechanical shock waves, such as those expected 
to be generated by a collision of two neutron stars, can 
cause repetitious superconducting phase transitions in 
nuclear matter. The relaxation of the proton conden- 
sate to its potential minimum can result in the non- 
thermal resonant production of 7-quanta in the MeV en- 
ergy range. 

There are some important consequences. The presence 
of energetic gamma quanta in nuclear matter, relatively 
transparent to photons thanks to electron degeneracy, 
can affect the equation of state. This, in turn, can affect 
the dynamics of the neutron star coalescence. 

In addition, the gamma rays trapped until the onset of 
the fireball but then released by the explosion, can also 
contribute to gamma-ray bursts. 
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